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Abstract 

The effect of the off-ramp (way out), on the traffic flow phase transition is in- 
vestigated using numerical simulations in the one-dimensional cellular automaton 
traffic flow model with open boundaries using parallel dynamics. When the off- 
ramp is located between two critical positions i c \ and i C 2 the current increases with 
the extracting rate /?o, for (3q < /3q c i, and exhibits a plateau (constant current) for 
0od < A) < Poc2 and decreases with @q for (3o > floc2- However, the density un- 
dergoes two successive first order transitions; from high density to plateau current 
phase at (3q = /?o c i ; and from average density to the low one at (3q = (3qc2- in 
the case of two off-ramps located respectively at i\ and 12, these transitions occur 
only when %i — i\. Phase diagrams in the (a, (3o),((3, (3q) and (ii,@o) planes are 
established. It is found that the transitions between free traffic (FT), Congested 
traffic(CT) and plateau current (PC) phases are of first order. The first order line 
transition in (ix, /?o)phase diagram terminates by an end point above which the 
transition disappears. 
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1 Introduction 



In the past few years, traffic problems have attracted the interest of a community 
of physicists [1-3]. Traffic flow, a kind of many-body systems of strongly interacting 
vehicles, shows various complex behaviors. Numerous empirical data of the highway 
traffic have been obtained, which demonstrate the existence of qualitatively distinct 
dynamic states [4-7]. In particular, three distinct dynamic phases are observed on 
highways: the free traffic flow, the traffic jam, and the synchronized traffic flow. It 
has been found out experimentally that the complexity in traffic flow is linked to 
diverse space-time transitions between the three basically different kinds of traffic 
[4]- 

To understand the behavior of traffic flow, various traffic models have been 
proposed and studied, including car-following models, gas-kinetic models, hydrody- 
namic models, and cellular automaton (CA) models [8-17], these later models have 
been applied to the various traffic problems for the automobile on the highway and 
the streets. First proposed by Nagel and Schreckenberg [9] and subsequently stud- 
ied by other authors using a variety of techniques. In the (CA) models, cares are 
treated as distinguishable particles and roads are expressed by discrete lattices. The 
whole system evolves in discrete time steps. All cars can move ahead simultaneously 
at each time step. These treatments for the behavior of cars make analyses of the 
complex systems with many degrees of freedom. With the help of these models, free 
flow and jams are well understood. On the other hand, the nature of synchronized 
traffic flow remains unclear despite various efforts [4,5,16,18]. 

Recent experimental investigation shows that in the majority of cases, synchro- 
nized traffic is observed localized near bottlenecks and thus it is believed that bot- 
tlenecks are important for the formation of synchronized traffic. The bottlenecks 
include on-ramps, off-ramps, lane closings, up-hill gradients, narrow road sections, 
etc. Among the various types of bottlenecks, the on-ramp is of particular interest to 
researchers and has been widely studied [16,18-21]. Popkov et al. [22] compared the 
experimental data from an on-ramp with simulations of a cellular automaton model. 
It is concluded that the dynamics due to the ramp can be described by a bottleneck. 
Diedrich et al [19] and Campari and Levi [20] simulated the on-ramp using the CA 
models. Traffic phenomena such as synchronized flow, the lane inversion, and phase 
separation were reproduced. 

The effects of the way(s) out (off-ramps) have been largely studied (empirically 
and numerically) in several recoveries. Our principal aim in these works is to answer 
to equation "where we can place the way(s) out and with which rate of absorption in 
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order to have a better fluidity? Using numerical simulations, we study the interaction 
between the traffic flow on the main road and way(s) out, i.e. the effect of the 
extracting rate (3q, in both one way out (one off-ramp) and two ways out (two off- 
ramps) in a road, on the dynamical jamming transition, density , current, and phase 
diagrams in the cellular automaton of one-dimensional traffic flow model with open 
boundaries using parallel update. It is found that the current exhibits a plateau for 
special positions of the way out. In spite of the simplicity of our model, it could have 
replicated the different behavior as well as phases fundamental of the traffic flow. 
CA models are conceptually simpler, and can be easily implemented on computers 
for numerical investigations. 

The paper is organized as follows; in the following section we define the model, 
the section 3 is reserved for results and discussions, the conclusion is given in section 
4. 

2 Model 

Motions of cars and interactions between cars are the microscopic processes in the 
traffic flow. One of the approaches to microscopic traffic processes is based on 
cellular automata. A CA model treats the motion of cars as hopping processes on 
one-dimensional lattices. The Wolfram's rule 184 CA is the simplest choice [23]. 
Nagel and Schreckenberg introduced their CA model (NS model) by extending the 
184 CA to consider the high velocity and the stochastic processes [9]. They showed 
that start-stop waves appear in the congested traffic region as observed in real 
freeway traffic. In this paper, we use the exclusion asymmetric models to serve as 
the basis model for the implementation of way(s) out. 

We consider a one-dimensional lattice of length L. Each site either occupied by 
one particle or is empty. A configuration of the system is characterized by binary 
variable Tj where Tj = (ji = 1) if site i is empty (full). To discuss the implementa- 
tions of way(s) out (off-ramp (s)), we assume that the main road is single lane and 
way(s) out connects the main road only on one lattice i\ in the case of one way out, 
and two lattices i\ and ii in the case of two ways out. During a time step At, each 
particle in the system has a probability p of jumping to the empty adjacent site 
on its right (and does not move otherwise), a vehicle can enter without constraint, 
with a probability a, in the first site being to the left side of the road if this site is 
empty. While, a vehicle being on the right in the last site can leave the road with 
a probability (3 and removed on the way(s) out (off-ramp (s)) with a probability /Jo- 
in our simulations the rule described above is updated in parallel (At = l.)[17], 
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i.e. during one update step the new particle positions do not influence the rest and 
only previous positions have to be taken into account. In order to compute the 
average of any parameter u (< u >), the values of u(t) obtained from 5 x 10 4 to 10 5 
time steps are averaged. Starting the simulations from random configurations, the 
system reaches a stationary state after a sufficiently large number of time steps. In 
all our simulation we averaged over 50-100 initial configurations. The advantage of 
parallel update, with respect to sub-lattice or sequential update is that all sites are 
equivalent, which should be the case of a realistic model with translational invari- 
ance. 

Thus, if the system has the configuration T\(t) ,T2(t) ,...,TL{t) at time t it will change 
at time t + At to the following, in the case of one off-ramp: 

For 1 < % < L and % ^ ii, n(t + At) = 1 with probability pi = Tj(t) + [Tj_i(t)(l — 
n(t)) - n(t)(l - r i+1 (t)]At and n(t + At) = with probability 1 - p { . 
Where At=l for parallel update. 

For i = 1, n(t+At) = 1 with probability p 1 = T 1 (t) + [a(l-T i (t))-T i (t)(l-T 2 (t)]At 
and T\(t + At) = with probability 1 — p\. 

For i = i\ and/or i = i 2 , Tj(i + At) = 1 with probability pi = Ti(t) + [rj_i(t)(l — 
n(t)) - /3 Ti(t)]At and n(t + At) = with probability 1 - p { . 

Where i\ and (3q denote the position(s) of off-ramp (s) and the extracting rate, re- 
spectively. 

Fori = L, T L (t+At) = 1 with probability p L = T L (t) + [T L -i(t)(l-T L (t))- pT L {t)]At 
and Ti(t + At) = with probability 1 — pl- 

3 Results and Discussion 

In the case of one way out ( one off-ramp) and in order to show the combinative 
effect of a, i\ and (3q on the traffic behavior, the variation of the average density 
as a function of {3q for several values of a is represented in Fig. la, we distinguish 
three regions when a is varied, a region in which we have a low value of density 
and does not depend upon (3q when a < j3, and when increasing a, high density 
and plateau current (PC) (the current undergoes a plateau as a function of (3q ) 
take place. We note that the MC phase appears between two critical values of (3q 
(A)ci an d A)c2)> and disappears completely when a becomes greater than a critical 
value contrary to the case of the one dimensional asymmetric exclusion model with 
open boundaries with parallel dynamics (At = 1.) [17] in which the MC phases 
disappears completely. However, the system exhibits three phases namely free traffic 
phase (FT), congested traffic (CT) and plateau current phase (PC); the transitions 
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between these phases are first order because of the discontinuity of the global mean 
density at the transition. The variation of the global mean density as a function of 
the extracting rate [3q for a = 0.4, (3 = 0.1, and for various positions of way out is 
given in Fig.lb. It is clear that the system exhibits three different phases, namely, 
FT, CT and MC. Indeed, there exist a critical position of the way out i c \ below 
which the density is constant for (3q < (3q C 2 and decreases discontinuously when 
increasing [3q for [3q > (3q C 2 • When the way out is located at a position greater than 
i C 2, the density decreases with (3$ for (3$ < [3q c \ and becomes constant for any value 
of /3o greater than I3q c i- While, when the way out is located at any position between 
i c i and i C 2 (i c i < i C 2), the density decreases, for (3q < (3od, with (3q and becomes 
constant for /?o c i < A) < A) C 2 and decreases with (3$ for (3q > (3q C 2- In the later 
case the density undergoes two successive first order transitions; from high density 
to intermediate one (plateau current phase) at /?o = (3q c \ and from intermediate 
density to the low one at (3q = (3qc2- In order to show the flow behavior under the 
effect of q and the positions of off-ramp, Fig.2a present the variation of the average 
flow as function of f3o for several values of a with (3 = 0.1 and i\ = 200, we show that 
the traffic flow increases with a, passes through a maximum at (3q = (3 max for high 
values of a and decreases for any value of (3q greater than (3 max . Indeed, for a certain 
values of a the current increase and becomes constant between two values of (3q and 
decreases for (3$ less than a critical values. Fig. 2b give for fixed values of a and 
(3, an indication of how much the average current change as ii is modified. When 
i\ < i c \ The current remain constant for (3q < (3q C 2 and decrease when increasing 
(3q. Moreover for (3q < [3q c \ the current increase and becomes constant for i\ > i C 2. 
While, when the way out is located at any position between i c \ and i C 2, the current 
undergoes three successive regions when (3$ is varied, a region in which the current 
increase for (3q < (3q c \ and remain constant between /3o c i and /3oc2 and decreases 
for (3q > (3q C 2- The system studied exhibits three phases: FT, CT and PC. The 
nature of the transition between these phases depends upon the values of a, (3, (3q 
and ii. Collecting the results illustrated in Fig.l and Fig. 2 we obtain the phase 
diagrams as shown in Fig. 3. The low and high density phases are separated by a 
first order transition line. Each of these phases undergoes discontinuous transition 
to the plateau current phase as show in Fig. 3a. To illustrate the effect of (3 for fixed 
value of a we give the phase diagram in the (3$) plan in Fig. 3b, we show that exist 
three phase similar in precedent phases and the transition between them is also the 
first order type. In the same way, one established another phase diagram but this 
time in the space (ii, (3q) as show in Fig. 3c. In the case of one way out, and for a fixed 
value of a and (3 the system exhibits three phases, namely, congested traffic, free 
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traffic and plateau current phases. The first order transition between the jamming 
phase and the moving one occurs at A) < 0.72 and the way out position lower than 
i& = 100. The same nature of transition occurs at A < 0.37 and the way out 
position greater than i& = 350. While, when the way out is located at any position 
between i c \ and i C 2, the density undergoes two successive first order transitions; 
from high density to intermediate one (plateau current phase) at A) = 0.37 and 
from intermediate density to the low one at A) = 0.72. For a larger values of a the 
plateau current disappears and the system exhibits a first order transition between 
high and low densities phases as Fig.3d and Fig.3e. In the case of Fig.3e the line 
transition terminates by an end critical point above which there is no distinction 
between the phases. In order to show the space-time evolution and the different 
region maps inside the chain. Moreover, the spatial vehicle density in the case of 
one off-ramp placed in the middle site of the chain depends on the values of a, (3 
and Pq. However for a = 0.4, (3 = 0.1 the system is in high (low) density phase for 
low (high) values of A; while for intermediate value of A) = 0.6 the coexistence of 
the free flow and jam takes place, as it is shown in the space-time evolution (Fig.4). 

In order to investigate the contribution that the extracting rate, Ab affects upon 
the density inside the chain. Fig. 5 gives for fixed values of a and (3, an indication 
of how much the average occupation of each position i (p(i) =< Tj >) changes as 
A is modified. The observed shifts may be interpreted as follows. For low values 
of A) (A) = 0.2), the average occupation is constant in the right side in the chain 
and will undergo a jump spontaneously in the position of the off-ramp and remain 
constant in the left one, for intermediate values of A (A) = 0.5), the value of average 
occupation remains constant for any positions less than the position of off-ramp and 
increase smoothly in the high value of average occupation then it becomes constant, 
for very large value of A) (A) = 0.8) p(i) remains unchanged for low position less than 
the position of off-ramp and decrease discontinuously in low value of occupation. 
We observe the same behavior in the case of two ways out with two jumps to the 
level of the position of the ways out. 

In the case of two ways out located symmetrically with regard to the center 
respectively at the positions i\ and 12, when the distance between two ways out is less 
than a critical value d c the traffic flow increases with Ab an d becomes constant at a 
critical value A)c' an d decreases at an other critical value A)c" ; increasing the distance 
between two ways out, the current increases with A), passes through a maximum at 
A = Pmax and decreases for any value of A greater than (3 max . However, the two 
successive transitions occur only when the distance 12 — i\ separating the two ways 
is smaller than a critical distance d c , and the system exhibits an inversion point at 
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/3o = /?0i = 0.4 for different values of distance between two ways out, as show in 
Fig.6a. Indeed, with increasing the distance 12 — h the current decreases for (3q > [3^ 
, while it increases for (3q < [3^. The inversion point disappears when the position of 
the first way out is located at the middle site; hence, the current increases with (3q 
and/or the distance ii — ii, passes through a maximum at (3q = (3 max and decreases 
for any value of (3q greater than (3 max (Fig.6b). The values of Pod, Poc2, id, V2 and 
d c depend on the injecting rate a, the extracting rate (3 and the position(s) of the 
way(s) out in the road. Moreover i c \ and i C 2, depend on the size of the system. 

4 Conclusion 

Using numerical simulation method, we have studied the effects of the extracting 
rate /?o, in both one way out (one off-ramp) and two ways out (two off-ramps) in a 
road, on the traffic flow phase transition, in the one-dimensional cellular automaton 
traffic flow model with open boundary conditions. We have shown that the behavior 
of density and current depends strongly on the value of (3q, the position of the way 
out from the entering and the distance between the ways out. In the case of one 
way out, and for a fixed value of a and (3 the system exhibits three phases, namely, 
moving phase, jamming phase and plateau current phase. The PC phase occurs 
only at special positions of the off-ramp (s) with special values of the extracting 
rate. The transitions between different phases are first order. Furthermore we have 
shown that the system exhibits an end point in the (off-ramp position, extracting 
rate) plane in the case of large value of the injecting rate. In the case of two off- 
ramps i\ and 22, the MC appears only when the distance between them is smaller 
than a critical value which depends on the absorbing, injecting and extracting rates. 
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Figure captions : 

Fig.l: The variation of average density p as a function of (3q, (a) %\ = 200 and 
for different values of a, (b) a = 0.4 and for different values of %\ , with /? = 0.1 and 
L = 400 

Fig.2: The variation of current J as a function of (3q, (a) %\ = 200 and for different 
values of a, (b) a = 0.4 and for different values of i\, with (3 = 0.1 and L = 400 

Fig.3 : Phase diagrams (a) (a, /So) plan for (3 = 0.1 and i\ = 200, (b)-(d) present 
phase diagrams in the {(3,(5q) plan for i\ = 200, a = 0.4 and a = 0.7 respectively, 
(c)-(e) phase diagrams in the (ii, /3o) plan for a = 0.4, (3 = 0.1, and a = 0.7, (3 = 0.2 
respectively, for L = 400 

Fig.4: space-Time evolution for alpha = 0.4, beta = 0.1, (3q = 0.6 and L = 400, 
in case of one way out. 

Fig.5: density profile in the case of one way out for a = 0.4, j3 = 0.1 and for 
several values of (3q 

Fig. 6 : The variation of the current versus (3q in the case of two ways out for 
different values of the distance 12 — h between the two ways, with a = 0.4, (3 = 0.1 
and L = 400, (a) the position of the two ways out change symmetrically with regard 
to the center of the chain, (b) the first way out is fixed in the middle site and the 
second change the position 
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